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ABSTRACT 

We present results of a 3D model of optical to 7-ray emission from the slot 
gap accelerator of a rotation-powered pulsar. Primary electrons accelerating to . 
high-altitudes in the unscreened electric field of the slot gap reach radiation- 
reaction limited Lorentz factors of ~ 2 x 10 7 , while electron-positron pairs from 
lower-altitude cascades flow along field lines interior to the slot gap. The cur- 
vature, synchrotron and inverse Compton radiation of both primary electrons 
and pairs produce a broad spectrum of emission from infra-red to GeV energies. 
Both primaries and pairs undergo cyclotron resonant absorption of radio pho- 
tons, allowing them to maintain significant pitch angles. Synchrotron radiation 
from pairs with a power-law energy spectrum from 7 = 10 2 — 10 5 , dominate the 
spectrum up to ~ 10 MeV. Synchrotron and curvature radiation of primaries 
dominates from 10 MeV up to a few GeV. We examine the energy-dependent 
pulse profiles and phase-resolved spectra for parameters of the Crab pulsar as a 
function of magnetic inclination a and viewing angle (, comparing to broad-band 
data. In most cases, the pulse profiles are dominated by caustics on trailing field 
lines. We also explore the relation of the high-energy and the radio profiles, as 
well as the possibility of caustic formation in the radio cone emission. We find 
that the Crab pulsar profiles and spectrum can be reasonably well reproduced 
by a model with a = 45° and ( ~ 100° or 80°. This model predicts that the slot 
gap emission below 200 MeV will exhibit correlations in time and phase with the 
radio emission. 

Subject headings: pulsars: general — radiation mechanisms: nonthermal — ac- 
celeration of particles — stars: neutron — X-rays: stars — gamma rays: theory 
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1. INTRODUCTION 


The study of emission from rotation-powered pulsars will soon undergo a major advance 
with observations by two new 7 -ray telescopes. The AGILE mission launched in April 
2007 and carries a pair-production telescope with sensitivity at 100 MeV comparable to 
the EGRET telescope on the Compton Observatory, but with much larger effective area 
at energies above 1 GeV. AGILE is expected to find 7 -ray pulsar counterparts to some of 
the unidentified EGRET sources. The Large Area Telescope (LAT) on the Gamma-Ray 
Large Area Space Telescope (GLAST), due to launch in early 2008, will have unprecedented 
sensitivity (30 times better than EGRET) and energy resolution for 7 -rays in the range of 
20 MeV to more than 300 GeV. GLAST is expected to significantly advance our knowledge 
of particle acceleration and radiation from rotation-powered pulsars, by discovering many 
more 7 -ray pulsars and measuring pulse profiles and phase-resolved spectra to much higher 
precision. Specifically, the LAT may increase the population of 7 -ray pulsars to several 
hundred and obtain flux measurements with error bars five-times smaller than those of 
EGRET. In order to make the best use of the power of LAT observations, one needs the 
most sophisticated and diverse models available. 

This paper adds to this effort by presenting, for the first time, results of a full model 
of high-altitude radiation from a pulsar slot gap and its relation to the geometry of the 
radio emission. Even 40 years after the discovery of pulsars, the origin of the particles 
that radiate high-energy photons as well as how and where they are accelerated is still not 
known. A number of models have been proposed and studies over the years, including 
polar cap models, putting the site near the neutron star at the magnetic poles (Arons & 
Scharlemann 1979, Daugherty & Harding 1982), and outer gap models, placing the site of 
acceleration and emission in the outermost magnetosphere near the light cylinder (Cheng et 
al. 1986, Romani 1996, Hirotani 2001). Such models have had varying degrees of success, 
with the polar cap models being able to explain very well the phase-resolved spectra of Vela 
-like pulsars (Daugherty & Harding 1996) and the outer gap models being able to explain 
more easily the double-peaked profiles of the Vela, Geminga and Crab-like 7 -ray pulsars. 
However, even multiwavelength observations have not been able to definitively rule out any 
of the proposed models. 

Studies of the geometry of polar cap accelerators (Arons k, Scharlemann 1979, Arons 
1983) discovered the possibility of a slot gap, a narrow bundle of field lines bordering the 
closed-field region where the parallel electric field is decreasing and accelerating particles 
cannot create pairs. The pair formation front across the polar cap thus forms a bowl shape 
as pairs are produced at higher and higher altitude approaching the last open field line. 
Although the existence of a slot gap has been known for some time, its potential as a 
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high-energy emission site has only recently been explored. Muslimov & Harding (2003) 
derived the electric potential and field in the slot gap, and modeled the acceleration and pair 
cascades at altitudes from the neutron star surface up to several stellar radii. They found 
that cascades develop over a range of 3 - 4 stellar radii along the slot gap inner boundary. 
These cascades have higher pair multiplicity (number of pairs per primary electron) than 
cascades nearer the magnetic axis because the primary electrons accelerate and produce pairs 
over a larger distance. The radiation from the slot gap cascades naturally produces a wider 
hollow cone of emission than the cascades near the surface. Muslimov & Harding (2004) 
extended the electric field solution in the slot gap to high altitudes and modeled only the 
curvature radiation of the primary particles. They found that such emission formed caustic 
patterns, like those of the “two-pole caustic” model proposed by Dyks & Rudak (2003) to 
naturally explain the double-peaked profiles of 7-ray pulsars. However, the spectrum of 
curvature radiation from radiation-reaction-limited particles alone is too hard to match that 
of observed 7-ray pulsars. 

The radiation model of Muslimov & Harding (2004) was not complete since it did not 
include all sources of primary particle radiation and all particles capable of contributing 
to high-altitude radiation. We have now included additional radiation components of the 
primary particles as well as radiation from electron-positron pairs, produced in polar cap 
and slot gap cascades at lower altitude, that continue to flow along open field lines with 
relativistic velocity. We investigate both the synchrotron radiation of both primaries and 
pairs, that undergo cyclotron resonant absorption of radio photons and significant increases 
in pitch angle, and non-resonant inverse- Compton scattering of radio photons. We also 
explore the geometry of the radio emission in relation to that of the high energy emission. 
In Section 2, we describe the geometry of the magnetic field and of the slot gap, Section 3 
discusses the model of acceleration in the slot gap and Section 4 details our treatment of 
the emission processes we model. We describe the simulation of radiation in the slot gap in 
Section 5 and present the results for the Crab pulsar in Section 6. Finally, in Section 7 we 
discuss the results and how the characteristic features of the model could be measured by 
GLAST and AGILE. 


2. 3D GEOMETRY OF THE MAGNETOSPHERE 

2.1. Retarded Vacuum Dipole and Open Volume Coordinates 

The modeling of the neutron star magnetospheric geometry is based on the calculation of 
Dyks, Harding & Rudak (2004, DHR04), who adopted the retarded vacuum dipole solution 
(e.g. Romani & Yadigaroglu 1995, Cheng et al. 2000) to account for the rotational distortions 
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and sweepback of the magnetic field lines. While the vacuum dipole solution does not 
accurately apply to a real pulsar magnetosphere that is filled with charges, and more realistic 
ideal MHD solutions now exist (Spitkovsky 2006, Timokhin 2006), the MHD solutions are 
numerical and time-consuming to simulate. The vacuum dipole solution exists in analytic 
form, is thus much easier to implement and we believe it exhibits the features like sweepback 
of field lines that approximate the MHD solutions. 

Our code first performs Kunge-Rutta integrations along the distorted field lines to de- 
termine the location of the open field footpoints at the neutron star surface that form the 
distorted polar cap rim. Details of this calculation are given in DHR04 and Dyks & Hard- 
ing (2004). We then define “open volume coordinates” (r ovc ,l ovc ) at the stellar surface to 
identify footpoints of field lines on which particles are traced. The “radial” coordinate r ovc 
ranges from 0 at the magnetic pole to 1 at the polar cap rim. Lines of constant r ovc thus 
form a set of distorted concentric rings that cover the polar cap surface (see Figure 2 of 
DHR04). The “azimuthal” coordinate l ovc measures the arclength along each deformed ring 
of fixed r OV c- The arclength is measured in the direction of increasing magnetic azimuth <j> pc , 
counterclockwise around the polar cap, starting from l ovc = 0 at 4> pc = 0, defined to be at 
the magnetic meridian (i.e. line between the magnetic and spin axes). One notable feature 
of the distorted polar cap rim is the appearance of a “notch” at </> pc ~ 30°, first discovered 
by Arendt & Eilek (1998) and most prominent between inclination angles of 40° and 60° 
(see Figure 2 of Dyks & Harding 2004). The notches are caused by field lines crossing near 
the light cylinder because of the sweepback, and produce bunching of fields lines at lower 
altitudes. This bunching of field lines will produce enhanced radiation at certain azimuth 
angles and noticeable features in the radiation phase plots, as will be discussed in §6. 

The open volume coordinates (‘ovc’) enable us to selectively model different particle 
distributions across the polar cap, such as electrons accelerating and radiating in the slot 
gap. To model a distribution of particles that is uniform in magnetic azimuth, we define a 
set of rings between r™|. n and r™** with equal spacing d ovc . The footpoints are then placed 
uniformly around each ring with = A azim 360 equal divisions. Since we trace emission 
from particles on the same number of magnetic field lines around each ring, and the rings 
have different circumferences, the contribution from different rings is weighted by Ling/Lim, 
where i ring is the ring length and / r i m is the polar cap rim length. 


2.2. Geometry of the Slot Gap 

The slot gap, a narrow region bordering the last open magnetic field line where particles 
cannot accelerate rapidly enough to produce pairs before the field drops, is an unavoidable 
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feature of polar cap space-charge limited flow (SCLF) models (Arons & Scharlemann 1979 
[AS79], Arons 1983). SCLF models assume that charges are freely emitted from the neutron 
star polar cap surface and flow out along open field lines. Since this charge flow is not 
sufficient to supply the Goldreich-Julian charge above the surface, an E\ \ exists and the 
charges are accelerated. Radiation from these charges forms electron-positron pairs in the 
strong magnetic field, which can screen the E\\ above a pair formation front (PFF) in a 
distance small compared to the acceleration distance, which is typically less than a stellar 
radius (AS79, Harding & Muslimov 2001). On field lines well inside the polar cap rim, E\\ is 
relatively strong and the PFF is very near the neutron star surface. These models assume a 
boundary condition that the accelerating electric field and potential vanish at the last open 
field line. Near the boundary, the electric field is decreasing and a larger distance is required 
for the electrons to accelerate to the Lorentz factor needed to radiate photons energetic 
enough to produce pairs. The pair front thus occurs at higher and higher altitudes near the 
boundary and curves upward, asymptotically approaching the last open field line (Harding & 
Muslimov 1998), forming a narrow slot gap. Since E\\ is unscreened in the slot gap, particles 
continue to accelerate and radiate to high altitude along the last open field lines. The width 
of the slot gap, A£sg, is a function of pulsar period and surface magnetic field, and can be 
expressed as a fraction in colatitude £ = 0/0 o of the polar cap opening angle, 0 0 ~ (flR/c ) 1/2 
(Muslimov & Harding 2003), 

« 2 P(AR 12 )- 4 / 7 / 45 3/7 . (1) 

where we will take A = 0.1. For the Crab parameters, rotation period P = 0.033, surface 
magnetic field Bn = .Bo/10 12 = 8 and moment of inertia I 45 = // 10 45 g-cm 2 = 4, the slot 
gap width A£ sg ~ 0.04. Here and in the rest of the paper, we will adopt these values, where 
Bn comes from the measured P and P, and A is a free parameter within the range 0.1 — 0.5 
that sets the width of the slot gap. Although the slot gap width itself is somewhat weakly 
dependent on these parameters, the slot gap potential and electric field are more strongly 
dependent on these parameter values (see §3). 

MH03 found that pair cascades on the inner edge of the slot gap occur at altitudes of 
3-4 stellar radii and have higher multiplicities M + ~ 10 4 — 10 5 than the polar cap cascades. 
Since the slot gap is very narrow for young pulsars having short periods and high fields, 
the corresponding solid angle of the gap emission Q SG oc 0 2 t]A(sg is quite small at low 
altitudes, 77 <C i]lc, where r/ = r/R is the dimensionless radius. However, this approximate 
expression does not describe the solid angle of radiation from the high-altitude slot gap, where 
relativistic effects strongly distort the emission pattern (see §6) into narrow caustics, and the 
solid angle must be computed by 3D mapping on the sky. So even though only a small fraction 
of the polar cap flux is accelerated in the slot gap, the radiated flux <f> S G = Lsg/^sg d 2 can 
be substantial. The particles in the slot gap can achieve very high Lorentz factors which at 
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altitudes of several stellar radii are limited by curvature radiation losses, to 7 sg — 3 — 4 x 10 7 
(Muslimov & Harding 2004). In this paper, we model the emission from the electrons in the 
slot gap as they continue to accelerate to high altitude above the region of pair formation. 


3. ACCELERATION OF PARTICLES IN THE SLOT GAP 


The solution for the electric potential and parallel electric field in the slot gap was 
presented by Muslimov & Harding (2003, MH03) in the low altitude limit and by Muslimov 
& Harding (2004, MH04) in the high altitude limit. Both solutions assume a space-charge 
limited flow of electrons from the neutron star surface, with the boundary condition p = p GJ 
at the stellar surface, where p is the true space charge density and p G j is the Goldreich-Julian 
charge density. The space charge density at the neutron star surface that includes relativistic 
frame-dragging, from Eqn (4) of MH03, is 


Po = 


QBci 


2 nccto 


[(1 — k) cosa + -9qH(1) since cos <f> pc }, 


( 2 ) 


where B 0 is the surface magnetic field strength in Gauss, Q, = 2 n/P and R are the neutron 
star rotation rate and radius, a is the pulsar obliquity, 0 pc is the magnetic azimuthal angle 
and 6 0 « [QR/cf( 1)] 1/2 is the polar cap half-angle. Here, k, = (r g / R)(I / 7 0 ) « 0.15 I 45 /RI is 
the general relativistic parameter from the frame-dragging effect, r g is the NS gravitational 
radius, R$ = R/ 10 6 cm and ao = (1 — r^/R) 1 ^ 2 . In addition, /( 1) ~ 1.4 and H( 1) ~ 0.8 are 
general relativistic parameters. Since the p(r) decreases faster than p G j(r) above the surface 
(for the explicit expressions see MH03), the deficit ( p — Pgj ) increases with distance, causing 
a parallel electric field, V • E\\ = Air(p — p G j) to develop. The treatment of the electric field 
in the slot gap by MH03 and MH04 differs from that of Arons & Scharlemann (1979) and 
Arons (1983) in that MH03 take into account the screening of the electric field by pairs on 
field lines interior to the slot gap (see §2.2). This forms fully conducting boundaries on both 
the inner and outer edges of the slot gap, causing a lower electric field in the slot gap. From 
Eqn (58) of MH04, the low-altitude solution for E\ \ is 


Ed, low —3(1 


UR \ 2 B 0 


c J /(l) 


K/ 1 

— cosa + -9oH(l)S(ri) sin a cos </> pc 
V ^ 


(3) 


where 6 (r/) (Muslimov k Harding 1997) varies between ~ 0.5 and 1, and v SG = is a 

parameter related to the width of the slot gap, A£ so , given by Eqn (1) above. The solution 
in Eqn (3) is valid for radii p < ?p, where ?/ c is a free parameter to be determined by matching 
smoothly to the high-altitucle solution. The high-altitude solution for valid for ?/ ~ ?/ c , 
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is given by Eqn (53) of MH04, 

3 /QR \ 3 


E i 


.high 


8 


Bi 

c > /(l) 


iV~{ 


, 1 /r 8 

1 + ( 5 3 

3 V Vc 


+ 2 


V_ 
Vic J 


COS Q+ 


—^o^7(l) since cos (f) pc 


(4) 


It is interesting that frame-dragging’s dominant effect on the accelerating field persists even at 
large distances from the neutron star surface, since the high-altitude SCLF solution depends 
on surface boundary conditions. Since a pure dipole field was assumed to derive the above 
expression for E\\, it is probably not accurate for very high altitudes approaching the light 
cylinder (r ~ 0.7/?lc), where the field lines will be distorted by effects of rotation and particle 
inertia. We combine the low and high-altitude solutions by use of the general expression 


E\\ a £j|,i ow exp [-(77 - l)/(i 7 c - 1)] + £j|,hi g h 


(5) 


and determine r] c to give a smooth transition between E|| >high and Ey ow . For the Crab pulsar, 
we find that t] c = 1.2 gives a smooth transition at most inclination angles. 


An estimate of the full potential drop of the slot gap can be obtained by s ummi ng that 
of the low altitude potential (Eqn [12] of MH03) from the surface up to r} c = 1.2 and that of 
the high altitude potential (Eqn [52] of MH04) at the r]i c : 

$Tot = *£foc) + *hgh (Vic) = Soflfoo 0.5(1 + k) cos a, (6) 

where $0 = (flR/c)B 0 R. In the case of the Crab pulsar, taking Bo = 8 x 10 12 G, Q = 190, 
v SQ = A£ 2 g /4 = 4x 10 -4 and K = 0.14, ~1.3x 10 13 eV. 

The particle flux from the polar cap and slot gap is 

A — cA sur , (7) 

e ' 

where 

^.ur = xROtiK T) 2 ~ (C?) 2 ] - xR9l£^ eo (8) 

is the surface area either of the slot gap or of the field line surface footpoints of the pairs, 

with p = p 0 for the primary electrons and p = A4 pairs (r 0 vc)po f° r the pairs. Again for the 
case of the Crab, hs G ~ 7 x 10 32 for the slot gap. 


4. EMISSION PROCESSES 

We model radiation over the entire spectrum from radio to 7 -ray wavelengths. The 
radio emission is necessarily treated phenomenologically, since no consensus physical model 
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has yet emerged. We therefore adopt a geometrical description of the radio emission beam 
as consisting of a core component and single conal component, described below. For the 
emission from optical through high-energy 7-ray, we simulate the radiation from primary 
electrons accelerating in the slot gap and also the radiation from secondary electron-positron 
pairs, that are produced in cascades near the neutron star surface and are flowing on field 
lines interior to the slot gap. 


4.1. Radio beam model 


We model the radio emission beam using an empirical core and cone model that has 
developed over the years through detailed study of pulse morphology and polarization char- 
acteristics. The average-pulse profiles are quite stable over long timescales and typically 
show a variety of shapes, ranging from a single peak to as many as five separate peaks. The 
emission is also highly polarized, and displays changes in polarization position angle across 
the profile that often matches the position angle swing expected for a sweep across the open 
field lines near the magnetic poles in the Rotating Vector Model (Radhakrishnan V Cooke 
1969). Rankin’s (Rankin 1993) study of pulse morphology concluded that pulsar radio emis- 
sion can be characterized as having a core beam centered on the magnetic axis and one or 
more hollow cone beams also centered on the magnetic axis surrounding the core. Although 
Rankin’s model assumes that emission fills the core and cone beams, other studies (Lyne & 
Manchester 1988) conclude that emission is patchy and only partially fills the core and cone 
beam patterns. 


The particular description we adopt is from Gonthier et al (2004) and is based on work 
of Arzoumanian et al. (Arzoumanian, Chernoff & Cordes 2002, ACC), who fit average-pulse 
profiles of a small collection of pulsars at 400 MHz to a core and single cone beam model 
based on the work of Rankin. The summed flux from the two components seen at angle 0 
to the magnetic field axis (modified by Gonthier, Van Guilder & Harding (2004) to include 
frequency dependence Sr) is 

S(0, Sr) = iwe-* 2/ ^e + F cone e- (e - §)2 ^ (9) 


where 


F,(sr) 


^-(l + I/,;)/ Sr \"i + l Li 

Sr \50MHz/ VtiD 2 


( 10 ) 


and the index i refers to the core or cone, v, is the spectral index of the total angle-integrated 
flux, Li is the luminosity of component i and d is the distance to the pulsar. The width of 
the Gaussian describing the core beam is 


Pc ore = 1.5°P-°- 5 


(11) 
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where P is the pulsar period in seconds. The annulus and width of the cone beam are 


9 = (1. — 2.63 5 w )p c , 

^wPconc 

where 6 W = 0.18 (Gonthier et al. 2006), and 


n — 1 94° r 0 - 5 P -0 - 5 
ycone — ±-a*± i kg r 


where 


0.07 


r K c ~ 40 


p0.3_-0.26 
r b GHz 


^10“ 15 ss _1 j 

is the radio emission altitude in units of stellar radius (Kijak & Gil 2003). 
The luminosities of the core and cone components are then 


Tcone 


Tradio j ^rad'io 

") l-'core 


1 + (r/r 0 )’ 


1 + (r* 0 /r) ’ 


where 


To 


G C one {y core + 1) 1 / £r 


Gcore (r'cone + l)r \s 0 


^core ^cone 


( 12 ) 

(13) 

(14) 

(15) 


(16) 


(17) 


where v coxe = —1.96 and u cone = —1.32, and 

L radio = 2.87 x 10 10 p- 13 p 0A mjy • kpc 2 • MHz 
where P is in units of lss -1 . The core-to-cone peak flux ratio is (Gonthier et al. 2006) 


(18) 


r = 


-{ 


OP pl.3_— 0.9 
zor f GHz> 


4 p - 18 £-°- 9 


GHz’ 

-0.9 

GHz’ 


P < 0.75 
P > 0.7 5 


(19) 


According to Eqn (15), the altitude of the conal radio emission is a weak function of the 
pulsar period but the emission occurs increasingly close to the light cylinder as the period 
decreases. For Crab-like periods, the conal emission occurs at altitudes of 10% — 20% of the 
light cylinder radius. 

We will use this radio emission beam model both to calculate the radio emission pro- 
file and to compute the intensity and angles of radio photons for the cyclotron absorp- 
tion/synchrotron emission component described in Section 4.3 below. To incorporate this 
radio emission geometry in the retarded dipole magnetic field that we are using to simu- 
late the high energy emission, we modulate the field lines of the ovc with the flux S{6, e R ) 
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given by Eqn (9). The differential flux radiated from a bundle of field lines centered at ovc 
coordinates (r ovc , / ovc ) is 

dS(»,e„) = S,(e,€„) S in9i vt «„r”“^<fe R (20) 

where iVj is the number of azimuthal divisions of each ring. The flux is assumed to be 
emitted at altitude 1.8 R for the core component and at altitude given by Eq (15) for the 
conal component. 


4.2. Curvature Radiation of Primary Electrons 


Primary electrons accelerating in the slot gap will produce curvature radiation up to 
7-ray energies. The curvature radiation losses will be balanced by the acceleration due to 
E\\ so that the particles reach steady-state Lorentz factors. The curvature radiation energy 
spectrum from a single electron with Lorentz factor 7 is 

N C r(£) = V3-'7k(—) (21) 

where e is the radiated photon energy in units of me 2 and 


t C r 



( 22 ) 


and the function k(x) is defined as 

roc 

k(x) = 2x K 5/3 (x')dx'. 
J 2x 


(23) 


The field line radius of curvature p c is not the standard one for a pure dipole field, but is 
determined in the ininertial observer’s frame by interpolating from tables computed from 
variations of the vacuum retarded field lines (see §2.1). The approximate form of the photon 
spectrum is a power law with an exponential cutoff at e cr , 


Ncr(s) 

de 


<*f 

(A me) 1 /' 3 \p c 


2/3 

•£“ 2/3 exp(-£/£ cr ). 


(24) 


4.3. Cyclotron Resonant Absorption and Synchrotron Radiation 

Relativistic particles moving parallel to magnetic field lines can gain perpendicular mo- 
mentum by undergoing cyclotron resonant absorption of low-energy photons. The cyclotron 


- 11 - 


resonant absorption of radio emission by relativistic particles in pulsar magnetospheres, fol- 
lowed by spontaneous synchrotron emission, was first proposed some years ago by Shklovsky 
(1970) as a mechanism for generating the optical radiation from the Crab pulsar. The process 
involves the absorption of photons at the cyclotron resonant frequency in the rest frame of 
the particle, resulting in an increase in the particle pitch angle. The particle then spon- 
taneously emits cyclotron or synchrotron radiation, depending on whether its momentum 
perdendicular to the magnetic field, in the frame in which the parallel momentum van- 
ishes, is non-relativistic or relativistic. Blandford & Scharlemann (1976) computed the cross 
section for cyclotron resonant absorption, but their application of the process to the Crab 
pulsar resulted in a re-radiated cyclotron radiation flux that was too small to explain the 
Crab optical emission. However, they assumed that the perpendicular momentum remained 
non-relativistic, in which case the applicable rate is that of cyclotron emission from the first 
excited Landau state, which is small relative to the rate from highly excited states. Lyubarski 
& Petrova (1998, LP98) performed a more detailed analysis of the distribution functions of 
particles undergoing synchrotron resonant absorption of radio photons and found that the 
particles can increase their pitch angles rapidly enough in the outer magnetosphere to attain 
relativistic perpendicular momentum. They found that the pitch angle excitation rate due to 
resonant absorption is much higher than the de-excitation rate due to synchrotron radiation 
until the particles reach very large pitch angles. The resulting synchrotron radiation ex- 
plain the optical and X-ray emission spectrum of the Crab and other young pulsars (Petrova 
2003). Harding et al. (2005, HVM05) showed that cyclotron resonant absorption of radio 
emission can work very efficiently for millisecond pulsars, especially for those pulsars where 
the accelerating electric field is unscreened. 

We apply the formulation of HUM05 to model the synchrotron radiation components 
from both primary electrons accelerating in the slot gap and from non-accelerating pairs 
on field lines just inside the slot gap. In the case of the primary electrons, the calculation 
closely parallels that of HUM05, who found that accelerating particles undergoing cyclotron 
resonant absorption will reach a steady-state where the synchrotron radiation losses are 
balanced by the acceleration gain. The resonant absorption condition is 

B' = 7£ 0 (1 - P/j.) (25) 

where 7 is the particle Lorentz factor, £0 is the energy of the radio photon in the lab frame 
(in units of me 2 ), (3 = (1 — 1 / 7 2 ) 1/2 , B' = B/B cr is the local magnetic field strength in units 
of the critical field B CI = 4.4 x 10 13 G, // = cos 6, and 0 is the angle in the lab frame between 
the photon direction and the particle momentum (to good approximation the direction of 
particle momentum is the same as the magnetic field direction). If this condition is met, 
the radio photon energy is at the local cyclotron energy in the particle rest frame. The 
resonant condition is not satisfied near the neutron star surface, where the magnetic field 
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is too strong, but may be achieved when the particle reaches the lower magnetic fields at 
high altitudes. When a particle is undergoing absorption initially in low Landau states, the 
rate of cyclotron emission is much lower than the rate of absorption. The Landau state and 
pitch angle of the particle will therefore increase continuously until an equilibrium is reached 
between gain in pitch angle through resonant absorption and the loss in pitch angle through 
synchrotron emission. Since this equilibrium is achieved at high Landau states, the emission 
is synchrotron rather than cyclotron. 

Lyubarski & Petrova 1998 identified two regimes of resonant absorption as the pitch 
angle of a particle increases. When ip <$; 9 (i.e. the pitch angle of the particle, ip is less 
than the incident angle of the radio photons, 9 ), the particle pitch angle increases but 
the total momentum stays roughly constant. When (9 — ip) «C 9, the pitch angle stays 
constant while the total momentum increases. Petrova (2002) has derived the solution for 
the distribution function of electrons undergoing resonant absorption of radio emission in a 
pulsar magnetosphere and the corresponding mean-square value of the pitch angle for both 
of these regimes. In the case where ip <C 9 (equation (2.17) of Petrova (2002)), the mean 
square of the pitch angle is 

(0 2 ) =4 R f a(jf)drf, (26) 

JVR 

where r? = r/R and 


27r 2 e 2 (l -0no)l o /eo7(l ~ @Po)Y 

( B' )' " > "- 


(27) 


Here Iq is the intensity of observed radio emission measured in erg ■ cm ~ 2 • s -1 • Hz~ l and v 
is the radio spectra index. Thus, for the perpendicular momentum change due to cyclotron 
resonant absorption we can write 


dp i 

dt 


abs 9 _ 

= 2 a 0 {ri) c— + — 
Pi P 



(28) 


where we used the relationship p _ L = p(0 2 ) 1,/2 . Thus, we assume that p±_ is proportional to 
the root mean- square value of the pitch angle. We also make the further approximation of 
computing the evolution of the root mean-square value of p± rather than the evolution of the 
particle distribution function. Since the primary and secondary electrons are continuously 
accelerating, 7 remains very high and Pxjp = sin ip <C 1. According to Petrova (2003), the 
width of the p± distribution is of order p± , so that the large variations in 7 and p in r along 
the field lines is much more important in the formation of the spectrum than the spread in 
the p± distribution. 
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Combining Eqn (27) and (28), the rate of resonant absorption can be written 


dpxY bs n Y Pi 7 (*Y Y bs 
dt ) Pi 7 2 - 1 \dt) ’ 


7 < 7 R 


(29) 


where 

D = 5.7 x 10 9 s _1 7r‘ / $o[mJy](l - 0fio) t (30) 

and we also can neglect the (d'y/dt) abs term in (29) since it is small compared to ( d'y/dt ) from 
acceleration and from curvature and synchrotron losses. In the above expression, we have 
assumed that Iq — $ 0 Cl rad d 2 /A, where T 0 is the measured radio flux (in mJy), d is the source 
distance, Q rad ~ A/r 2 is the radio emission solid angle, with A and r the cross-sectional area 
and radius at the absorption radius. Also, from the resonant condition (see equation [25]), 
7 < 7m Jr is defined as 


7 R = 


B' 


= 2.8 x 10 5 - 


B* 


£o(l — 0IM}) ^0,GHz(l — Pno) 

The resonant terms will switch on only when the resonant condition is satisfied. 


(31) 


In the case where 9 — ip <§; 9, we assume that the pitch angle remains at the value 
ip — 9/2 and the total mean particle momentum evolves as 


P !-(£) 


(3 - i/)aiY /ia ~ v) 

b x 9 2 


(Petrova 2003), where 


, , 47T 2 e 2 J'?/ 0 f eo9 Y 


(32) 


(33) 


and b\ = 2e 2 B l2 /3h 2 c. 

To evaluate the local intensity of radio photons 4> 0 and the incident absorption angle /i 0 
needed in Eqns (30) and (33) we use the radio core/cone beam model described in Section 
4.1 above. To compute To, we represent the emission at each point of the core and cone 
beam in ovc coordinates as a “beamlet” with opening angle /j, bm centered on the tangent to 
the field line and sum the contributions from all the separate beamlets. The beamlets that 
are modulated by the core beam are located at 1.8 R and the beamlets that are modulated 
by the cone beam are located at r KG (see Eqn (15)). To evaluate n Q for each beamlet that 
will encounter a particle at position ( x p ,y p ,z p ), we must take into account the aberration of 
the radio photons at their emission point (x bm ;y brn ,z bm ) as well as the light travel time from 
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the radio photon emission point to the particle and the rotation of the field line (and of the 
particle on that field fine). Because of the time delays and the rotation of the magnetic field, 
not all beamlets will encounter a particle at a given position in the outer magnetosphere. 

At each step along a field line (see §5), the particle radiates an instantaneous synchrotron 
spectrum, given by (Tademaru 1973) 

o2/ 3 

Nsr(s) = r— jr OLfB' sint/>£ - 2 / 3 £“p / 3 exp(— £/£ sr ), (34) 

1 \ 3 ' 

where sini/> = p±/p, p 2 = 7 2 - 1 and s SR = (3/2)7 2 B' sin^ is the synchrotron critical 
frequency. It is important to point out that the main radiation power produced by the 
particles undergoing resonant absorption of radio photons comes mostly from the parallel 
energy of the relativistic particles, not from the power of the observed radio emission (which 
is relatively much smaller) . The absorption of radio photons increases the pitch angles of the 
particles which already have very high 7 . For high 7 , the pitch angle remains approximately 
constant while the particle radiates synchrotron radiation, because the radiated photons 
are emitted nearly perpendicular to the magnetic field in the frame of circular motion, and 
therefore nearly along (within angle ~ I/ 7 ) the particle momentum in the lab frame. The 
particle then experiences recoil in a direction opposite to its motion, decreasing 7 x and 7 || 
in proportion to each other, allowing the parallel component of energy to be tapped. 


4.4. Non-resonant inverse Compton 


The primary electrons will also scatter radio photons that are not in the cyclotron reso- 
nance in their rest frames. We treat this component as a non-resonant Compton scattering, 
which will be in the Thomson limit since 7£ 0 < 1. The spectrum of scattered photons can 
be written 


dNjcs 

d£ s dfx s 


J d *J ie C imAe ' ^ 


(35) 


where £ and e s are the incident and scattered photon energies in the lab frame, p and p s 
are the cosines of the incident and scattered photon angles in the lab frame and n P h(s,/i ) 
is the number density of incident radio photons. The primes denote the corresponding 
quantities in the electron rest frame and are related to those in the lab frame by the Lorentz 
transformations 


J 


A*' 


7 £ (1 - 3 11) 
/i — 3 


(36) 

(37) 
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with same type of expression for the scattered quantities. The differential cross section 
for scattering of an electron in a magnetic field, averaged over photon polarization, in the 
Thomson limit, is written (e.g. Dermer 1990) 


da'(£',n') 3 j „ , 

= - Gjd{e a —e a ) 




Ms 2 )(1 — M /2 ) + ^(1 + a4 2 )(1 + A^ 2 ) 


tr 


+ 


i r 


L(u+1) 2 (u-l) 2 J 

(38) 

where u = s'/B 1 . Since the primary electrons have high Lorentz factors, the incident photons 
will be beamed into a narrow cone with O' ~ l/q so that n' ~ 1. The radio photons will also 
lie at energies near or above the cyclotron resonance in the electron rest frame (see §4.3), 
and since we are treating only the non resonant part of the scattering here, can approximate 
the cross section as : 




(39) 


For the incident radio photon distribution local to the primary electrons, we use the form 


n ph {£ , A-0 = n R d{e - £r)Q(/j,) (40) 

where n R is the local radio photon density, we approximate the energy distribution as a 5- 
function at energy £ R . We estimate the radio photon density as 

T ra di 0 / A1\ 

Ur £ 0 A{r)c 

where comes L radio from Eqn (18) and A(r) ~ irr 2 [(6+w e ) 2 — (0— w e ) 2 ], assuming the electron 
interact mostly with photons from the conal component. For the angular distribution we 
adopt 

(1 — n 2 ) cos 2 01 


9(a0 = 


1 

( 1 - P R /i ) 3 


1 - 


7|(1 -Pr/J .) 2 J 


(42) 


which is that expected for emission from relativistic particles with Lorentz factor 7# (Jackson 
1965). Equation (42) is meant to describe the distribution of coherent radio emission from 
pairs with Lorentz factors in the range 7« ~ 10 2 . In order to calculate the spectrum of 
scattered photons, we change variables from /j, to e' in Eqn (35) by means of the Lorentz 
transform in Eqn (36). The 0 integration can be done easily. We then make use of the 5- 
functions to perform the e and s' integrations and perform the integration over the scattered 
photon angle n s numerically. 


5. SIMULATION OF MULTI WAVELENGTH SLOT GAP RADIATION 

In order to simulate the radiation on the sky for a non-rotating distant observer, we have 
adopted two coordinate systems. The acceleration and emission takes place in the corotating 


16 


frame (CF). The frame of the distant, non-rotating observer we will refer to as the inertial 
observer frame (IOF) where the radiation is observed. We start the simulation by computing 
the distorted polar cap rim at the neutron star surface in the CF. The foot-point of the field 
lines on which the primary electrons are accelerated are positioned along concentric rings 
in the ovc and the (xo,yo,zo) coordinates along each field line are computed. The primary 
electrons are confined to the slot gap, which we assume is bounded by the ovc coordinate 
rings r™' c n = 0.95 and r™** = 1.0 with equal spacing d ovc = 0.004 giving 11 rings. Each ring 
is divided into 180 equal divisions. Each electron starts at the stellar surface on a trajectory 
that advances along the field line in increments determined by either its rate of acceleration, 
radiation loss rate, perpendicular momentum gain or distance gain. The electron trajectory 
extends to a radius r max = 0.8 R\ c , since the field fine structure of the retarded vacuum 
dipole solution is not accurate very near the light cylinder. However, the resulting emission 
pattern and profiles do not depend very sensitively on emission patterns in the outermost 
parts of the magnetosphere, since the bulk of emission originates at lower altitudes (see § 6 ). 
At each step along the field line in the CF, the equations of motion for the electron (Eqns 
[43] and [44]) are solved to give the new Lorentz factor and perdendicular momentum. The 
spectrum of curvature, synchrotron and ICS photons radiated in that step are computed and 
the radiation direction is assumed to be tangent to the local field direction in the CF. The 
emission direction is then transformed to the IOF (aberration), appropriate time delays are 
added and the radiation from curvature, synchrotron and ICS are accumulated in separate 
arrays P(e,C, </>), where ( and <f> are the viewing angle and phase with respect to the pulsar 
rotation axis for a distant, non-rotating observer. 

The electron-positron pairs are also treated in a similar way as the primary electrons, but 
with the acceleration gain rate set to zero, since they are created above the pair formation 
front which forms the inside boundary of the slot gap where the parallel electric field is 
screened. The pairs are assumed to flow along open field lines bounded by the ovc coordinate 
rings r™ n = 0.25 and = 0.99 with equal spacing d ovc = 0.014 giving 54 rings. On rings 
r OV c < 0.25, the field line radius of curvature increases near the magnetic axis and the pair 
multiplicity rapidly decreases (Daugherty & Harding 1996), so we ignore pair radiation on 
these very inner rings. We ran our code on a parallel processor, where each ring was assigned 
to a different processor. 


5.1. Particle dynamics and radiation 


The Lorentz factors, 7 , and perpendicular momentum, p ± (in units of me), of each 
particle will evolve along the field lines according to its equation of motion, which may be 
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written 


dr) eE\\ 2e 4 2 2e 2 7 4 /'cky'N 

P± ~ ~3rt + V dt ) 


(43) 


dpx = _3c 2 e 4 p2 p\ f dpjj) \ 
dt 2r 1 3 m 3 c 5 7 \ dt ) 


(44) 


The terms of the right hand side of equation (43) are acceleration, synchrotron losses, cur- 
vature radiation losses and cyclotron/synchrotron absorption. The terms of the right hand 
side of equation (44) are adiabatic changes along the dipole field line, synchrotron losses 
and cyclotron/synchrotron resonant absorption. A derivation of the above equations (minus 
the CR loss and resonant absorption terms) is given in the Appendix of Harding, Usov & 
Muslimov (2005). The ICS losses may be neglected for the primary particles and may also 
be neglected for the pairs since the acceleration and synchrotron loss rates are much larger. 


By substituting equation (28) into the right hand sides of equations (43) and (44), we 
get 

^ = A i E ll -B 1 Blpl-C n i * (45) 

Q—Mn-'K-KBiA (46) 

where A x = 1.76 x HFs" 1 , B t = 1.93 xlOV 1 , Ci = 5.6 x lO^s" 1 , A 2 = 4.5 x 10 4 s" 1 , E\\ is 
in e.s.u. and Bg = B/10 8 * * * G. Since Petrova (2002) has assumed that p and 7 are constant to 
compute the change in pitch angle due to resonant absorption, we have neglected the change 
in 7 due to absorption in equation (45). Both E\\ and Bg are functions of rj. For the primary 
electrons, E\\ is computed using Eqns (3), (4) and (5). For the pairs, we assume E\\ = 0. 


6. RESULTS FOR THE CRAB PULSAR 

The Crab, with a rotation period P = 33 ms, and period derivative P = 4.22 x 10 ~ xz s/s, 

implying a surface magnetic field of B 0 = 8 x 10 12 G, is the youngest pulsar having detected 

pulsed 7-rays. The 7-ray profile shows two sharp peaks with a phase separation of about 0.4 

or 144°. Very similar profiles appear at all wavelengths, including radio, where an additional 
pre-cursor peak leading the first or main peak is present at frequencies below about 450 

MHz. We model the Crab spectrum and pulse profile using our model of emission from 

primary electrons accelerating in the slot gap and non-accelerating electron-positrons pairs 
flowing along field lines inside the slot gap (i.e. field lines at smaller colatitude). Polar cap 
pair cascade models (Daugherty & Harding 1996; Muslimov & Harding 2003) find that pair 
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multiplicity (number of pairs per primary particle) for the Crab is in the range 10 4 — 10 5 , 
where the lower number applies to colatitudes nearer the magnetic axis while the latter 
number applies to cascades extending to higher altitudes of 3 — 4 stellar radii, initiated 
by electrons accelerating in the slot gap. In simulating the radiation from pairs, we trace 
particles along all open field lines with an assumed broken power law spectrum of energies 



7_ min < 7 P < 7 


br 


•P 

w br 


Tp <7 p < % 


(47) 


where the normalization constants are set to preserve the total pair multiplicity A4 pa irs(fovc) 
in each ring. We explore two geometries for the radio emission beam: the ‘standard’ beam 
model described in §4.1 with the core emission at a fixed altitude and cone emission at 
a single frequency-dependent altitude, and a non-standard model with core emission at the 
same fixed altitude but the cone emission is extended along the last open field line to produce 
caustic peaks. 


Figure 1 plots several examples of the evolution of particle Lorentz factor and per- 
pendicular momentum that are solutions to the equations of motion (43) and (44) along 
magnetic field lines at leading and trailing edges of the open volume, for the extended ra- 
dio cone emission case. Also plotted are the acceleration gain rate, d'y at:c /dt, curvature loss 
rate, d^/^/dt, synchrotron loss rate, d'jsjdt , cyclotron absorption rate, d / y a , bs /dt, and critical 
synchrotron energy t SR . In the top plots, showing these quantities for primary electrons, 
the electric field acceleration falls off rapidly from its initially high value near the neutron 
star surface but maintains a nearly constant lower value out to high altitudes. The electron 
Lorentz factor is limited by curvature radiation losses at fairly low altitude, with the gain 
rate from acceleration balancing the curvature loss rate along nearly its entire path (thus 
the two curves lie on top of each other). The absorption begins at the lowest altitude of 
conal radio emission, r m ; n ~ 0.2 Ri c at 1 GHz, producing a sudden increase in p±, which then 
increases steadily with altitude. As a result, the synchrotron losses turn on and reach a level 
nearly that of curvature losses on the leading field lines, but somewhat lower on trailing 
field lines. The fluctuations in absorption rate are due to the electron encountering radio 
emission from different beamlets (see §4.3) and so reflects the numerical resolution of these 
sub-elements of the radio cone beam. At higher altitude, the radio cone emission turns off, 
the electrons on trailing field lines no longer encounter radio photons and p± begins to drop 
along with synchrotron losses. The synchrotron radiation reaches energies of several hundred 
MeV. The bottom plots show the evolution of these quantities, except for acceleration gain 
(which is zero) and curvature losses (which are negligible), for an electron or positron of 
Lorentz factor 7 p = 10 5 near the upper end of the pair spectrum (see below). The evolution 
of pi and cfyabs/df show similar behavior to that of the primary electrons, except that the 
pairs ‘see’ the radio photons out to higher altitude since they are on more interior field lines. 
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The synchrotron emission of these highest-energy pairs reaches energies of a few MeV, while 
that of the lowest energy pairs with = 10 2 peaks at optical frequency. 

We adopt a model where the slot gap high-altitude emission best reproduces the Crab 
phase-averaged spectrum. This model has M = 1.4 M 0 , R = 14 km, I = 4 x 10 45 gcm 2 and 
A = 0.1 for the parameter values determining the primary electron acceleration. For the 
parameters of the pairs spectrum, we assume 7 ™ in = 10 2 , 7 ™ ax = 2 x 10 5 , ^ = 5 x 10 3 , 
<5i = 2.0, 82 = 2.8, which matches the pair distribution derived for the Crab parameters in 
Fig. 7 of Daugherty & Harding (1982). The pair multiplicity was varied in rings over the 
polar cap, such that 

{ 4 x 10 3 0.25 < r ovc < 0.5 

2 x 10 4 0.5 < r ov c < 0.9 (48) 

4 x 10 5 0.9 < r ovc < 0.99 

The multiplicity of cascades on the field lines near the magnetic axis is lower than that for 
cascades near the slot gap. Since the E\\ near the slot gap has a lower magnitude, primary 
electrons accelerate and produce extended cascades over much longer distances. The primary 
electrons also produce most of the cascades during their radiation- reaction limited phase, 
further increasing the multiplicity of pairs near the slot gap. 


6.1. Standard radio beam 

Using the output arrays P(s, £, <f>) from the simulations, we can display plots of radiation 
intensity on the sky, or phaseplots, in different energy bands. Figure 2 shows phaseplots of 
high-energy emission assuming the standard radio emission model for three energy bands: 
1 — 20 keV, 0.1 — 10 MeV, and > 100 MeV. The radiation distribution shown in the high- 
energy phaseplots exhibits caustics, extended bright lines of emission from particles on the 
trailing field lines from each pole. The near cancellation of phase shifts due to retardation, 
aberration and field line curvature cause emission at a wide range of altitudes to arrive in 
phase, while emission on the leading edge of the open volume is spread-out in phase (Morini 
1983, Dyks & Rudak 2003). Pulse profiles for these energy bands are obtained by displaying 
the intensity as a function of phase at a particular viewing angle £. At viewing angles that 
cut across caustics, the profiles show two peaks with phase separation less than 180°. We also 
show the phaseplot at 400 MHz of the radio emission, which is dominated by the cone beam. 
The cone beam is shifted earlier in phase relative to the core beam due to the difference 
in aberration and retardation of the higher-altitude cone emission. We find that inclination 
angles in the range a = 40° — 55° give the best combination of spectrum and profile to match 
the Crab profile and spectra. For large inclination angles, the parallel electric field reverses 
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direction on some field lines (see MH03), and for small inclination angles, radiation from the 
low-altitude pair cascades dominates the emission and produces double-peaked profiles for 
small viewing angles (MH03). For inclination angle a = 45°, viewing angles in the range 
£ = 78° — 82° and £ = 98° — 102° produce profiles having two peaks with phase separation 
near 0.4. The phase-averaged emission spectra are very similar for the different angles in 
this range. The profiles for energy bands 1 — 20 keV, 0.1 — 10 MeV look identical because 
both are due entirely to pair synchrotron radiation whose geometry is not energy dependent. 
The profile for energies > 100 MeV are significantly different, with the first peak now larger 
than the second peak, and the phases somewhat shifted from those at lower energy. While 
the high-energy profiles can reasonably reproduce the observed Crab profiles, the standard 
model radio profile does not reproduce that observed. The high energy phaseplots and 
profiles are in units of phs _1 / s t er and ph s^/ster/N^, respectively, where N# = 180 are 
the number of phase bins. The phase-averaged flux for an observer at viewing angle ( is 
then the sum over the emission in the profile divided by source distance squared. For the 
viewing angle £ = 100° shown in Figures 2 and distance of 2 kpc, the phase-averaged flux is 
($(> 100MeV)) = 4.3 x 10~ 6 phs -1 cm -1 and ($(0.1 — 10MeV)) = O.OSphs^cm -1 , which 
are in good agreement with observed values. 

In Figure 3, we show the our model phase-averaged spectrum for the same case a = 
45°, C = 100° shown in Figure 1. Four components are visible and three components make 
significant contributions to the total spectrum. Synchrotron radiation from pairs contributes 
at the lowest energies, from infra-red and optical through hard X-rays, turning over at around 
20 MeV. The range of this component reflects the pair spectrum, with the low energy turnover 
determined by 7 ™ in , the high-energy turnover set by 7™ ax , the break dividing lower and 
upper slopes determined by 7 ^ r , <5i and 82 - Synchrotron radiation from primary electrons 
contributes at the mid-range of 20 — 300 MeV and has a smaller energy range, reflecting the 
smaller energy range of the primaries. Finally, curvature radiation from primaries contributes 
at the highest energies, from ~ 100 MeV to the turnover at around 5 GeV that is determined 
by the parameters of E \\ . The component due to non-resonant ICS from primaries and from 
pairs (which appears below the scale of the plot) make negligible contributions. The non- 
resonant Compton scattering contribution is much lower than that of resonant absorption 
because the cyclotron absorption cross section is orders of magnitude higher, being a first 
order process, than the Thomson cross section for non-resonant ICS, being a second-order 
process. Furthermore, the particles undergo absorptions of many photons at the resonant 
before radiating the most significant synchrotron emission, reaching high Landau levels. 
Thus the ratio of resonant absorption to ICS is even larger than the simple ratio of the cross 
sections. Curvature radiation from pairs is also negligible since their Lorentz factors are 
much lower than that of the primary electrons. Our model spectrum plotted in Figure 3 has 
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not been adjusted arbitrarily to fit the data. Rather, we have adjusted the model parameters 
within a reasonable range. The pair spectrum can be tuned to match the optical-to-X-ray 
spectrum very well by adjusting its energy range and shape, as well as the variation of 
multiplicity across the polar cap. However, the values that best fit the data are not far from 
those that come out of pair cascade calculations. The model spectrum does not match the 
high-energy spectrum quite as well, with the peak of the curvature spectrum that gives a 
high enough cutoff energy being somewhat above the data points, although the EGRET data 
points have large errors. The EGRET sensitivity above 1 GeV has recently been re-examined 
by Stecker et al. (2007), and they concluded that the Crab phase-averaged spectral points 
(and indeed those of all EGRET sources) should be systematically lowered by as much as a 
factor of 2 above 1 GeV. We also plot these corrected points in Figure 3 and they provide an 
improved match to our model spectrum. The multicomponent nature of the model spectrum 
produces dips around 20 and 200 MeV. The observed spectrum seems to show the dip at 
20 MeV but a dip at 200 MeV is not clear, given the large errors. Measurements with the 
GLAST LAT should define such details much better. 

The phase-resolved spectra are shown in Figure 4, for the phase ranges adopted by 
Kuiper et al. (2001) to study the Crab pulsed emission. Our model spectra are obtained by 
summing the emission for the equivalent phase intervals (our Peak 1 occurs at phase 0.41 and 
Kuiper et al.’s Peak 1 occurs at phase 0.0) in the phaseplots at the chosen viewing angle. 
The flux levels of the resulting model spectra were then plotted with the data from the 
corresponding phase intervals without any renormalization. Generally, the phase-resolved 
spectra do not fit the data as well as the phase-averaged spectrum. This may be due to 
the fact that the low-energy profile peaks (due to pair synchrotron emission) are not exactly 
phase-aligned with the peaks of the high-energy profile (due mostly to primary curvature 
emission). However, the data and model spectra match fairly well for the PI, P2 and bridge 
intervals. The primary synchrotron component makes a larger contribution in the Peaks and 
the Bridge, which is expected since the primary electrons radiate only along a narrow set of 
field lines. 


6.2. Extended radio cone beam 

The phase alignment of the first and second peaks in the profile across frequency bands 
from radio to high-energy q-ray is a hallmark of the Crab pulsar and difficult to explain 
unless the emission originates from the same location in the magnetosphere. As we have 
seen in the last section, the caustics that form in emission along the trailing field lines at the 
edge of the open volume naturally produce phase-aligned peaks, even if the emission does not 
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originate at the same location along the field lines. A requirement for caustic peak formation 
is emission that is extended along these field lines above altitudes of a few tenths of R [c . One 
possibility to explain the phase alignment of radio and high-energy peaks in the Crab profile 
is that the radio cone emission, which is already expected in the standard model to come 
from altitudes above 0.1 Ri c , is extended by several tenths of Ri c in altitude. To explore 
this possibility, we modify the standard cone beam model by spreading the emission from 
r m in = ^kg to fjnax = Nr kg and decreasing its width w e , preserving the same luminosity. 
To implement this extended radio beam in the cyclotron resonant absorption/synchrotron 
emission component, we divide the beamlets of incident radio photons among a number of 
discrete altitudes between r min and r max . 

The resulting phaseplots and profiles for the extended radio emission model are shown 
in Figure 5. The radio emission indeed forms caustics similar to those seen in the high-energy 
emission phaseplots. To produce radio caustics we need to extend the radio conal emission 
from r m i n = tkc to r max = 57 "kg and take its width as 0.2iy e . The resulting radio profile 
at the same ( = 100° displays two peaks separated by 0.4 in phase that are at the same 
phases as the high-energy peaks. The high-energy phaseplots and profile calculations use the 
same parameters for the primary electrons and pairs as the model described in the previous 
section, the only difference being the radio beam geometry. There are only minor differences 
in the high-energy profiles, with a higher and more evenly distributed emission level in the 
bridge region and a less extended trailing wing on Peak 2 above 100 MeV. 

The resulting high-energy phase-averaged spectra, shown in Figure 6, are also very 
similar to those of the standard radio beam case in Figure 3. The level of the pair synchrotron 
spectrum is slightly lower but the shape is the same. The shape of the primary synchrotron 
spectrum is somewhat different, but the peak is still around 100 MeV. The primary curvature 
radiation is of course unchanged since it does not depend of the radio emission. There are 
somewhat larger differences in the phase-resolved spectra, shown in Figure 7, since they 
are more sensitive to the details of the radio emission geometry. The pair synchrotron 
components are lower in the leading wings, LW1 and LW2, but higher in the trailing wings, 
TW1 and TW2, and about the same in the peaks, PI and P2, and bridge intervals. Similarly, 
the primary synchrotron components are lower in the leading wings and higher in the trailing 
wings, with the exception of TW2, but higher in the Bridge. 


7. DISCUSSION 

A model of radiation from young, rotation-powered pulsars has been investigated, as- 
suming that the emission we observe occurs at high altitudes in the magnetosphere. We in- 
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clude two distributions of particles radiating along open field lines at high altitude: primary 
electrons accelerating in the slot gap from the stellar surface to near the light cylinder, and 
non- accelerating electron-positron pairs on field lines bordering and interior to the slot gap. 
Simulating spectral formation in 3D, we are able to study pulse profiles in different energy 
bands as well as phase-resolved spectra. We choose to first model the Crab pulsar, since it is 
bright in all wavebands from optical to high-energy 7-rays, with good phase-resolved spectral 
measurements. In our model, three different components dominate the Crab phase-averaged 
spectrum, curvature and synchrotron radiation of primary electrons and synchrotron radia- 
tion of pairs. The synchrotron radiation results from cyclotron resonant absorption of radio 
photons in the low magnetic field at high altitude. Cyclotron resonant absorption of ra- 
dio photons using standard core and cone emission models can reproduce high-energy pulse 
profiles but not the observed radio profiles or the phase alignment of high-energy and radio 
peaks. Extension of conal radio emission along the last open field lines results in caustic 
formation that can both produce a double-peaked radio profile and the phase-alignment with 
high-energy peaks. However, we find that the peaks in the model profiles at energies below 
100 MeV are not exactly phase-aligned with the peaks in the profile above 100 MeV. Also, 
our model phase-averaged spectrum does not reproduce exactly the observed spectrum above 
100 MeV, because the curvature radiation component appears high relative to the that of 
the primary synchrotron radiation. The observed spectrum in the peak intervals is matched 
better by the model spectra where the curvature radiation component is relatively lower. We 
predict that, due to the multicomponent nature of the model spectra, transitions between 
different components may be observed as dips in the peak and phase-averaged spectrum 
around 20 and 200 MeV which may be verified by GLAST LAT observations. 

Since absorption of radio photons is critical to the production of the synchrotron emis- 
sion components, this model predicts that the radio and high-energy emission in the syn- 
chrotron dominated parts of the spectrum, should be correlated in time as well as in phase. 
In our phase-averaged model spectrum of the Crab, the synchrotron radiationn dominates 
below about 200 MeV, so we predict that correlated variability should be seen in optical, 
X-ray and 7 rays below 200 MeV. Although Lundgren et al. (1995) detected no observable 
correlation between radio giant pulses and EGRET flux, Shearer et al. (2003) found small 
but significant (3 %) enhancement of optical flux during Crab giant radio pulses. Lorn men et 
al. (2007) found a correlation between the arrival time of pulses from Vela with X-ray profile 
shape as measured by RXTE. Such a correlation suggests a strong connection between radio 
and X-ray emisson. All of these results are consistent with our prediction if the EGRET 
photons that Lundgren et al. (1995) considered were dominated by the curvature radiation 
component (~ 200 MeV), which should not exhibit any correlations with radio emission. 

Although we have found from purely geometrical consideration that radio emission ex- 
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tended along the outer edge of the open field lines can produce a double-peaked profile whose 
peaks are in phase with those of the high-energy profiles, we must ask why the radio emission 
geometry of the Crab pulsar should be so different from that of other pulsars, including ap- 
parently the other known 7-ray pulsars. But if one looks at known X-ray pulsars, the radio 
and X-ray profiles of some millisecond pulsars such as PSR B1957+20 and PSR B1821-24 
also have phase-aligned peaks. In fact, the standard core and cone radio geometry is based 
primarily on emission morphology of pulsars having much longer periods in the range 0.5 - 
1 s. It is possible that the radio emission geometry in faster pulsars could be quite different. 
One could speculate that the slot gaps, that are predicted to exist only in the young and/or 
fast pulsars, sustain a beam of high energy primary particles near the last open field line 
that can interact with a flow of lower energy pairs, or an opposite flow of particles in a return 
current, to produce instabilities that lead to an extended radio component. These pulsars 
are also the ones that have giant radio pulses, a phenomenon that may be related to the 
non-standard emission geometry (Hankins & Eilek 2007). 

Takata & Chang (2007) and Tang et al (2007) recently presented results of a modified 
outer gap model for energy dependent pulse profiles and phase-resolved spectra of the Crab 
pulsar. In that work the classic outer gap, which extended from the null charge surface to the 
light cylinder along the last open field line, is assume to be extended below the null chanrge 
surface, so that the LW1 and TW2 parts of the pulse can be produced by emission from 
the second magnetic pole. The geometry of the modified outer gap model and that of the 
slot gap model presented in this paper are now quite similar. All the different parts of the 
profile come from the same regions of the magnetosphere in both models, with the exception 
of the first peak, which in the outer gap primarily comes from near the light cylinder. In 
the slot gap model, the first peak comes mostly from the caustic of the pole opposite to 
the one producing the second peak, although our > 100 MeV profile also has substantial 
components from both poles. The radiation mechanisms of the outer gap are very different 
from those of our slot gap model, with the major difference being the link to radio emission 
through resonant absorption of the radio emission in the slot gap model. While the optical 
to hard X-ray component of the Crab spectrum is produced by synchrotron radiation of 
pairs in both models, resonant absorption is responsisble for the pair pitch angles in the slot 
gap model but in the outer gap the pair pitch angles come from cascades outside the gap. 
Also, the 10 MeV to 500 MeV spectrum comes from inverse Compton emission in Tang et 
al (2007), while in our model this component comes from synchrotron radiation of primaries 
undergoing resonant absorption. The highest energy part of the Crab spectrum, above a 
few hundred MeV, comes from primary curvature radiation in both models, although this 
component originates only from very high altitude in the outer gap model. 

Although we have not presented results on the polarization here, this will be explored 
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in future work. The expected properties of polarization in a two-pole caustic geometry have 
been investigated by Dyks et al. (2004). They found that the position angle and polarization 
degree as a function of phase produce a reasonable match to the optical polarization of the 
Crab pulsar (Kanbach et al. 2005). Since the radiation modeled in this paper includes 
emission on interior field lines as well as in the slot gap, we expect that the polarization 
properties could be somewhat different and perhaps also energy dependent. 

Other Crab-like pulsars, such as PSR B0540-60 and PSR B 1509-58, have very similar 
X-ray to low-energy 7 -ray spectra to that of the Crab, but notably no EGRET detections. 
However, their single peaked high-energy profiles are quite different from the Crab profile, 
suggesting a different viewing angle. Both of these pulsars are at larger distances, so that the 
EGRET upper limits are above the level of emission they would have with an exact match 
of the Crab spectrum. A GLAST LAT detection of the high-energy emission component(s) 
would thus be able to probe the geometry of the slot gap or of other radiation models. 

These pulsars and also the Vela-like pulsars will be investigated in future work. Vela is 
predicted to have a slot gap so the spectrum above 100 MeV should be similar to that of 
the Crab. But since the pair multiplicity should be smaller than for the Crab, one expects a 
lower pair synchrotron component which may also peak at a different energy. The relation 
of the radio to high-energy peaks in Vela-like pulsars is quite different, with a single radio 
peak leading double-peaks at high energy. One possibility is that the radio cone is at too 
low an altitude to form caustic peaks (Harding 2005). 

In this work we have assumed steady-state acceleration in a space-charge limited flow 
model. However, since the time-dependence of such accelerators has not been investigated, 
this may not be correct. In fact, the magnitude of the current in SCLF models has not 
been determined self-consistently but is assumed to be constant at the Goldreich-Julian 
value across the polar cap. This current distribution is not consistent with that found in 
global force-free models of the pulsar magnetosphere (Spitkovsky 2006, Timokhin 2006), 
that require a current that departs from the Goldreich-Julian value in the outer parts of the 
polar cap. Matching the current distribution required by the global models with the average 
boundary condition on the charge flow of the SCLF accelerator implies either temporal or 
spatial fluctuations in the current and voltage. A lower primary charge flow in the slot gap 
or a smaller slot gap width may in fact produce a better match of the model spectrum to 
the observed spectrum of the Crab. Investigation of the temporal stability and boundary 
conditions of acceleration models will ultimately be needed in order to accurately model and 
understand pulsar high-energy emission. 
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Fig. 1 — Evolution of the dynamics and radiation of particles as a function of distance (in 
units if light cylinder radius) along leading (0 pc = 7r/2) and trailing (0 pc = -ir/2) magnetic 
field lines. Quantities plotted are particle Lorentz factor, 7 , perpendicular momentum, p ± 
(in units of me), acceleration gain rate, d^Jdt (s -1 ), curvature loss rate, dj CI /dt (s' 1 ), 
synclirotron loss rate, d%Jdt (s -1 ), cyclotron absorption rate, d^bs/dt (.s _1 ), and critical 
synchrotron energy s SR (in units of me 2 ). 
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Fig. 2. — Model pulse profiles and intensity maps (observer angle £ vs. rotation phase 
<j>) in different frequency ranges as labeled for the case of standard radio beam model and 
inclination angle of a = 45°. Profiles are shown for observer angle £ = 100°. The scale of 
the high energy intensity maps is in units of phs _1 /ster and the scale of the profiles is in 
units of ph s — 1 /ster/N^, where N# — 180 are the number of phase bins. The scale of the 
radio map is in mJy — kpc 2 /ster. 
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Fig. 3. — Model spectrum of phase-averaged total pulsed emission (heavy solid line) that is 
the sum of emission components from curvature (light solid line), synchrotron (dashed-dot 
line) and inverse Compton (dashed-dot-dot line) radiation of primary electrons in the slot 
gap and synchrotron radiation from pairs (dashed line) inside the slot gap. Data points are 
from Kuiper et al. (2001) [http://www.sron.nl/divisions/hea/kuiper/data.html]. The open 
squares are corrected EGRET values above 1 GeV from Stecker et al. (2007). 
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Fig. 4. — Model spectrum of total pulsed emission in different phase intervals (heavy solid 
line) that is the sum of emission components from curvature (light solid line), synchrotron 
(dashed-clot line) and inverse Compton (dashed-dashecl-dot line) radiation of primary elec- 
trons in the slot gap and synchrotron radiation from pairs (dashed line) inside the slot gap. 
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Fig. 5.— Model pulse profiles and intensity maps (observer angle C vs. rotation phase <p) 
in different frequency ranges as labeled for the case of an extended radio cone model and 
inclination angle of a — 45°. Profiles are shown for observer angle ( = 100°. 
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g. 7.— Same as Figure 4, but for the case of an extended radio cone model. 







